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We construct an effective Hamiltonian for electrons in the fractional quantum Hall regime for 
GaAs and graphene that takes into account Landau level mixing (for both GaAs and graphene) 
and sub-band mixing (for GaAs, due to the non-zero- width of the quantum well). This mixing 
has the important qualitative effect of breaking particle-hole symmetry as well as renormalizing the 
strength of the inter-particle interactions. Both effects could have important consequences for the 
prospect that the fractional quantum Hall effect at v = 5/2 is described by states that support 
non-Abelian excitations such as the Moore-Read Pfaffian or anti-Pfaffian states. For GaAs, Landau 
level and sub-band mixing break particle-hole symmetry in all Landau levels and sub-band mixing, 
due to finite-thickness, causes additional short-distance softening of the Coulomb interaction, further 
renormalizing the Hamiltonian. We find that in graphene, Landau-level mixing only breaks particle- 
hole symmetry outside of the lowest Landau level (JV 7^ 0). Landau level mixing is likely to be 
especially important in graphene since the Landau-level mixing parameter is independent of the 
external magnetic field and is of order one. Our realistic Hamiltonians will serve as starting points 
for future numerical studies. 

PACS numbers: 71.10.Pm, 71.10.+a, 73.43.Cd 



I. INTRODUCTION AND MOTIVATION 

One of the outstanding experimental challenges in 
physics is to determine whether the fractional quantum 
Hall effect supports non-Abelian anyon excitations. An 
affirmative answer would constitute the discovery of a 
new type of, so far unobserved, particle and could pave 
the way for topologically-protected quantum information 
processings^. 

The most promising and notable system for which 
the search for non-Abelian anyons is taking place is in 
the v — 5/2 fractional quantum Hall effecpHSI The 
search is motivated by two conjectures. The first is that 
the Moor e- Read Pfaffian stated and the anti-Pfaffian 
stat d 16 * 1 ^ are representatives of universality classes that 
have non-Abelian anyon excitations. This conjecture has 
been recently shown to be true^S (see also Refs. 1151 - 
[50]) . The second conjecture is that the experimentally 
observed state of matter responsible for the v = 5/2 frac- 
tional quantum Hall effect (FQHE) is in one of these two 
universality classes. Several numerical studies^^U have 
provided evidence supporting this conjecture by showing 
that the ground states of simplified model Hamiltonians 
are in these two universality classes (Moore- Read Pfaffian 
or anti-Pfaffian) . Even though these model Ha miltonians 
are physically reasonable, a number of studies 32 34 35 38 
highlight the sensitivity of numerical result a 31 " 35 * 37 '' 39 * to 
the parameters of these 'toy model' Hamiltonians (and 
to some degree, the system size). Experiments have 
found that relatively minor modifications to the system 
such as transposition of the half-filled Landau level from 
v = 5/2 to v = 9/2,11/2,... or to v = 1 /2,3/2 lead 
to metallic states that are anisotropic ! 40 * 41 ! or isotropic 



and Fermi liquid-like 4 ^. Similar anisotropic or isotropic 
metallic phases are also observe d at v = 5/2 upon ap- 
plication of, respectively, a smalP^Ell or large^ in-plane 
magnetic field (compared in magnitude to the perpendic- 
ular magnetic field). Thus, there are at least four differ- 
ent phases (and, perhaps, many more) which can occur 
at v = 5/2, depending on the details of the Hamiltonian. 
Thus, we are encouraged to ask basic questions which 
remain unanswered: is the ground state of a more real- 
istic model Hamiltonian for v = 5/2 in one of these two 
non-Abelian universality classes? What is the quantum 
phase diagram of a realistic Hamiltonian for v = 5/21 

The sensitivity of the v = 5/2 system to small changes 
in the Hamiltonian is a cause for concern since the canon- 
ical simplified model ignores both Landau-level mixing 
and the non-zero-width of the quasi-2D electron sys- 
tem. (The exceptions are Refs. I37|38l which include 
Landau-level mixing and Refs. 3 1 .'Hi which include non- 
zero-width.) Landau level mixing and non-zero-width 
can be neglected if the Landau level mixing parameter 
k = (e 2 /e£o)/huj c <C 1 and the non-zero-width of the 
quantum- well d/£o <C 1, respectively. Here, e 2 /e£o is 
the scale of the Co ulomb e nergy, uj c is the cyclotron fre- 
quency, and £0 — yjTic/eB is the magnetic length, so the 
dimensionless parameter k is given by k w 2.52/ -J 'B\T\. 
Therefore, for experiments at magnetic fields in the range 
IT— 10T, k is in the range 0.8 — 2.5. Similarly, for most 
GaAs samples d/£o « 2 — 3. It is not obvious that either 
one of these parameters can be considered small. Given 
that the ground state at v = 5/2 depends sensitively on 
the precise Hamiltonian and that the Hamiltonians stud- 
ied in all previous numerical results neglect potentially 
important effects, one may question their connection to 
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experiments- 6 

In fact, recently the FQHE at v = 5/2 has been stud- 
ied experimentally in the regime of strong Landau level 
mixing (ac > 2.5) and interesting non-linear behavior of 
the FQHE energy gap has been observed as a function 
of k (or density) 48 . Thus, in order to make meaningful 
experimental predictions for the FQHE in GaAs (or gen- 
erally any two-dimensional electron gas), it is necessary 
to study realistic Hamiltonians which include the effects 
of both Landau level mixing and non-zero-width. This 
is true not only at v = 5/2, but throughout the N = 1 
Landau level, where Landau level mixing is a factor. 

The scenario described above for GaAs is poten- 
tially an even more pressing issue when one consid- 
ers the FQHE in graphene. Due to the linear dis- 
persion of electrons in graphene, the cyclotron energy 
is sgn(N)y/2\N\hv F /£ , compared to tiu}(N + 1/2) for 
GaAs. (In GaAs, N = 0, 1, 2,. . . whereas in graphene 
— oo < N < oo). This means that the Landau level mix- 
ing parameter in graphene k — (e 2 / 'e£o) / '{-\f2Tiv p j ' 1$) — 
e 2 / {e\/2Tiv p) is, interestingly, independent of magnetic 
field strength. In addition, the spacing between Landau 
levels varies as l/\/2N for large iV-the Landau levels get 
closer together in energy the higher up, or lower down 
below N — 0, one goes. Consequently, the effective Lan- 
dau level mixing parameter for the A^ th Landau level is 
oc k\/N, which increases with increasing Landau level. 

Experimentally, the FQHE has been observed in 
graphene in suspended samplea^HH] (i. e ., free standing 
graphene) and on a boron nitride substrate^!. In sus- 
pended graphene, the Landau level mixing parameter is 
approximately k ~ 2.2 whereas on a boron nitride sub- 
strate it is much lower due to the reduction in dielectric 
constant: k w 0.44. In both cases k is obviously not 
small so it would appear that there is no experimental 
observation of the FQHE in graphene for which one can 
ignore Landau level mixing. Many theoretical works have 
considered the FQHE in graphene^HUl^ pointing out the 
similarities and differences between the FQHE in GaAs 
and graphene, but none have taken into account the ef- 
fect of Landau level mixing. 

Therefore, in this paper we construct and study a re- 
alistic Hamiltonian for GaAs in the lowest two Landau 
levels (N — and N = 1; no FQHE has been observed 
in N > 2) that takes into account Landau level and sub- 
band mixing. In addition, we construct a realistic Hamil- 
tonian for graphene in the N — 0,±1,±2 Landau levels 
which includes Landau level mixing effects which have, 
so far, been totally ignored. Although we are emphasiz- 
ing the experimental systems of GaAs and graphene in 
our calculations, they are appropriate for any fermionic 
system with either parabolic (GaAs) or linear (graphene) 
bands up to the specific experimental parameters. Once 
these effective Hamiltonians are characterized, they can 
be used to study a variety of problems in which sub-band 
and/or Landau level mixing effects may play a prominent 
role. 

Our starting point is the systematic treatment of 



Landau-level mixing formulated in Ref. 1531 The effec- 
tive Hamiltonian for electrons in the N th Landau level 
can be derived by integrating out all other Landau levels 
via an expansion in powers of A«(or k). This was done^ 
for GaAs for zero width, but the same procedure allows 
us to integrate out higher sub-bands in the same way 
as higher Landau levels and/or consider graphene, as we 
discuss below. 

Crucially, we find an error in the normal-ordering of the 
effective Hamiltonian of the previous work 59 . Correcting 
this error, we find that the renormalization of the two- 
body interaction is significantly modified and we examine 
the implications. 

We emphasize that, for GaAs, we include the non- 
zero width of the 2D layer, as a result of which the 
effective interaction between electrons is "softened" at 
short distances since the single-particle wavefunctions are 
smeared out over a length scale d in the direction perpen- 
dicular to the plane. This effect was found to stabili ze 
the Moore- Read Pfaffian and anti-Pfaffian states^^&l m 
the N = 1 Landau level - unlike in the lowest Landau 
level, where it weakens the FQHE. A further effect of 
non-zero-width is that electron-electron interactions can 
cause mixing with the higher quantum well sub-bands 
corresponding to motion in the direction perpendicular 
to the two-dimensional plane, in a manner analogous to 
Landau-level mixing. (Strictly speaking, a system with 
only sub-band mixing could also cause a breaking of 
particle-hole symmetry. However, Landau-level mixing 
is generally a stronger effect than sub-band mixing.) We 
take these two facets of non-zero width - the softening of 
the 'bare' Coulomb interaction and mixing with higher 
sub-bands - into account in a systematic manner. This 
has not, to the best of our knowledge, been previously 
done. Of course, for graphene, the effective width of the 
2D graphene sheet is negligible, so there are no sub-bands 
that need to be considered. 

Strictly speaking, our Hamiltonian can only be justi- 
fiably called 'realistic' for small values of k (or k). For 
GaAs, small k corresponds to higher values of the mag- 
netic field whereas for graphene, small k corresponds to 
a substrate with high dielectric constant. However, we 
find that the coefficients in the expansion in powers of k 
are small in GaAs (even with the correction mentioned 
in the previous paragraph) and consequently, our Hamil- 
tonian may even be realistic when k is not small. For 
graphene, the coefficients in the N = Landau level are 
small, similar to those in GaAs, but for N ^ 0, this 
is no longer the case. Thus, these expansions may be 
valid even in the regime where most experimental obser- 
vations of the FQHE at v = 5/2 in GaAs have taken place 
(0.7 < k < 1.8), as well as where the FQHE has been ex- 
perimentally observed in graphene^HSll (q.44 < k < 2.2). 
Our calculations might also help explain some of the pe- 
culiarities of the graphene FQHE observations, such as 
the fact that the FQHE has been observed in only the 
lowest Landau level, see Refs. HM5"2"1 This point will be 
discussed further below. 



3 



Previous results for GaAs: The Landau- level mixing 
Hamiltonian of Ref. [59] for GaAs was studied for zero 
width at v = 5/2 by Wojs et alW^ using exact diag- 
onalization in the spherical geometry, and it was con- 
cluded that, over nearly the entire range of Landau level 
mixing of experimental interest, < k < 3, the over- 
lap between the ground state and the Moore-Read Pfaf- 
fian state was larger than either the overlap between the 
ground state and the composite fermion Fermi-liquid-like 
wavefunctioipSl or the overlap between the ground state 
and the anti-Pfaffian wavefunction. However, the com- 
putation of Wojs et a/P^ did not take finite-thickness of 
the 2D system into account; it is not clear that these 
overlaps at different values of the magnetic flux track the 
(extrapolated) ground state energies, which is what de- 
termines the true ground state; and, more importantly, 
the Hamiltonian did not contain the corrected normal- 
ordering of the three-body term. Rezayi and SimorP^ 
did a similar study using the torus geometry but they 
did not use the Hamiltonian of Ref. [SS]and instead simu- 
lated Landau level mixing by diagonalizin g in a truncated 
Hilbert space. In contrast to Wojs et aZp3, they found 
the anti-Pfaffian to have a higher overlap with the exact 
ground state. Even though both studies used wavefunc- 
tion overlap to measure to which universality class the 
ground state belonged, it is difficult to directly compare 
their contrasting results since they used different Hamil- 
tonians and different geometries. Of course, it is possible 
that the previous error in the normal ordering of the 3- 
body term is the origin of these apparently contradictory 
results. 



Previous results for graphene: None of the theoreti- 
cal studies for the FQHE in graphene mentioned above 
included the effects of Landau level mixing. 



The plan of this paper is as follows. In Section |TT] we 
review the effective action description given previously 
in Ref. 59, To this we add the additional complication 
of a non-zero-width quantum well. This is not as simple 
as just augmenting the Coulomb interaction because it 
becomes necessary to integrate out all higher quantum 
well sub-bands. We then derive the effective Hamiltonian 
which follows from this effective action, paying particular 
attention to the operator-ordering of the Hamiltonian. 
In Section Hill we discuss the effective Hamiltonian which 
includes Landau-level and sub-band mixing to first order 
in the Landau level mixing parameter k. In Section IV 
we consider the effective Hamiltonian for graphene and 
in Secti on [V] we discuss the pseudopotentials. Finally, in 
Section | VI| we present our conclusions. 



II. EFFECTIVE HAMILTONIAN INCLUDING 
LANDAU LEVEL AND SUB-BAND MIXING 
APPLICABLE TO GaAs 

A. Effective Action 

We first begin by discussing the calculation for a sys- 
tem appropriate for GaAs, just to introduce the formal- 
ism, and then in Sec. |IV| we discuss the calculation for 
a system appropriate to graphene. Following Ref. 1591 we 
begin with the action for electrons in a magnetic field: 



S = 



/dcu 
2~ X! c snma (a;)(iw - E sn + n)c snma (uj) 
snm 

x C S2n2m2j g(w2)c Sl n imi a( l ^'l) 

X 27T(5(W4 + UJ3 — Ld2 — . (1) 



Here, 



and c s 



are Grassman variables where 



a =t, i are the spin indices; s = 0,1,... is the sub- 
band index; n — 0, 1, . . . is the Landau level index; 
and m — 0, 1, . . . labels orbital states within a Lan- 
dau level (771 — n is the angular momentum). E sn — 
huj c (n+l/2)+e s is the cyclotron energy (hcu c — heB/mc) 
plus the sub-band energy e s , and \x is the chemical po- 
tential. Note that the Zeeman energy has been set 
to zero in the above action since it is much smaller 
than the other energy scales in the problem. = 
(s4n4?7i4, 3377,37713! V\siri2mi, s\n\m\) is the interaction 
matrix element, where V is the electron-electron inter- 
action operator for the bare two-body Coulomb interac- 
tion for electrons in an infinite quantum well of width d, 
i.e. the usual non-zero- width augmented interaction in- 
vestigated previouslj'^lElISIHsS] Following the standard 
procedur e! 67 1 68 1 , we define G a b(k) as 



G a b(k) 




\k 2 /2) 



(2) 



for integers a > b (note that if a < b we switch k — 
k x + iky with k — k x — ik y ) with L b a (x) a generalized 
Laguerre polynomial. We find that the interaction matrix 
element is 



43.21 



cPk, 



Vs 4 s 3 s 2 s 1 (k)e G n4ri2 ( k)G n3ni (k) 



xC m4m2 ( k)G m37ni (k*) , 



(3) 



where the 4 in V^^i, for example, is short for all the 
internal degrees of freedom of the 4-th particle (34, 77,4, 
and 777,4) and 



V.. 



e 2 1 

2-7re/o k 



dzidz 2 (t>s i { z i) ( t>s 3 (zi) 

51 (^K fc|Z1 ^ 21 , (4) 
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where 4> s (z) is the wavefunction for the z-dependence tain: 
of single-particle wavefunctions in the s th sub-band. 
Throughout this work on GaAs we will consider the elec- 
trons to be confined to an infinite square well of width 
d. Thus, <j) s (z) = y / 2/dsm((s + l)nz/d) with z e [0,d] 
and e s (d) = (s + 1 ) 2 it 2 Ti 2 / (2m z d 2 ) where m z is the ef- 
fective electron mass in the quantum well. The total 
single particle energy is E sn = huj c (n + 1/2) + e s (d). 
If we assume that m z equals the electron band mass 
m and write d in units of the magnetic length £q — 
^/hc/eB w 25 nm/y/BlT], then E sn = hu c [(n + 1/2) + 
(■K 2 /2d 2 )(s + l) 2 ]. For typical quantum well widths and 
magnetic fields (d/£o ~ 2 — 4) the Landau-level spacing 
and sub-band spacing are comparable. Of course, the 
real potential is not an infinite square well, however, the 
infinite squ are w ell captures the key feature of the real 
potentiaP^EiimMl the ability of electrons to avoid each 
other by moving in the z-direction, as reflected both by 
the softening of the short-distance part of the Coulomb 
interaction and by the possibility of virtual transitions to 
higher sub-bands (which have nodes at certain values of 
z). 



1 f TT Aw% ( 2 ) - f \- i \ 

J 4=0 7 '' 

~ 3! /iT ^ U e%32lC0N7n 6 a(u6)c N7n 5 l3{u5) 

X Co Nm 4 7 ( w 4 ) CQJVro 3 7 i^V c 0iVm 2 /3 \ u 2) CONmi a (<^1 ) 
X27T 5{u>6 + + W4 — W3 — W2 — ^l) • (5) 

As a result of the higher Landau levels which are inte- 
grated out, the bare two-body interaction is renormalized 



,(2) 



'43 



;21 = ^43;21 + K^43;21 



(6) 



(3) 

and a three-body interaction Ug54-32i ^ s generated. At 
higher-orders in k, four-body, five-body, etc., interac- 
tions are generated, but we will restrict ourselves here 
to lowest-order in k, so we only need to consider two- 
and three-body interactions. 
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We first assume the ./V th Landau level to be partially 
occupied by electrons and integrate out all higher/lower 
Landau levels and sub-bands perturbatively in the 
Coulomb interaction. We thereby generate an effec- 
tive Hamiltonian for the ./V th Landau level which in- 
corporates the effects of virtual transitions to the other 
Landau levels and sub-bands. This effective Hamilto- 
nian is computed perturbatively in powers of K, where 
K = (e 2 /elo)/(hu] c ) is the Landau level mixing parame- 
ter discussed above (k « 2.52 x 1T/- V /S[T]). Thus, k 
can be varied from very small to very large by changing 
the magnetic field (and the density, in order to keep the 
filling fraction fixed) . This is is an important distinction 
between GaAs and graphcnc, as we discuss in more detail 
below. 



At tree level, the effective action is simply Eq. ([TJ) 
with the sums over Landau levels restricted to rij = N 
and the sums over sub-bands restricted to s, = 0. Most 
numerical studies of the quantum Hall effect use the ef- 
fective Hamiltonian associated with this effective action, 
in which the other Landau levels are integrated out to 
zcroth order in the Coulomb interaction or, equivalently, 
to zeroth order in k. However, if we integrate out the 
higher Landau levels to second-order, which amounts to 
computing the effective action to first-order in k, we ob- 




FIG. 1: The bare interaction Feynman diagram. The 
Coulomb interaction (augmented or not) cannot switch the 
spin so the first term on the right hand side has 8 aa S 1313 and 
the second term has 8 aP ' 8^ a . The labels 1, 2, 3, and 4 cor- 
respond to 1 = {si, m, mi}, etc.. 



The calculation will be presented using Feynman dia- 
grams and the single Coulomb vertex in the diagramatics 
is V& a = 7 4 3,2i<5 aQ 'SP? - V 3 4^i6 a ^'d^' where a, /?, 
a', and j3' label spin indices, see Fig. [T] 

We use the nomenclature of Shankar° 9 (as did Ref . I55|) . 
(see Fig. [2]) for the three diagrams which renormalize the 
two-body interaction as 



6u%t. n = ZS + ZS' + BCS 



(7) 



where ZS/ZS' stands for "zero-sound" and BCS stands 
for "Bardeen-Cooper-Schrieffer" due to their similarity to 
the corresponding diagrams in Fermi liquid theory. The 
corresponding expressions are: 



ZS EE 



00 dw, 

-oo 2 ^ 



2_ Z_^i v x3,x'l r 4aj',2a; 



X27T(5(W3 + UJ X — LJl 



)G X 'G : 



(8) 
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FIG. 2: (a) The ZS Feynman diagram: V^^ a v£^f 7 . (b) 
The ZS' Feynman diagram: V>? '?' a V?'?' ' f> ~ 1 . (c) The BCS 
Feynman diagram: Vf^'f 1 ' V^'f" ■ 



ZS' = - 



duj x 
~2tT 



v x4,x'l v 3x',2x 



, 27T ^ 

x2ttS(uj x + LU4 — w\ — uj x ')G x >G, 



(9) 



and 
BCS 



_i r°° du>x 

2 J-oo 27T 7-oo 27T ^ 

x2tt5(uj x > + lo x - oj\ - u) 2 )G x >G. 



v xx',21 v 43,xx' 



(10) 



where the electron propagator is G x — l/(iui x — E So , n:c ) 



with E s 



\i. We assume the chemical poten- 



tial to be the energy of the partially filled Landau level, 
/i = hoj c (N + 1/2), and the above S is shorthand for the 
sum 



E- E E' E' E 



(ii) 



Note the primes on the sums over sub-bands (s x and s x i) 
and Landau levels (n x and n x /) indicates that we do not 
simultaneously include the conditions (n x = N, s x = 0) 
or (n x i — N, s x > — 0)-we do not want to integrate out 
the N th Landau level of the lowest sub-band. 

We now turn to the evaluation of these expressions. 
Consider the ZS term, for example, and integrate over 
uj x i using the (5-function: 



ZS 



/OO J 
r, x3,x'l v 4x',2x 

-oo Z7T 



1 



i(u) 3 - wi + uj x ) - E Sm , n ^, ioj x 



Eg 



(12) 



and we can do the remaining integral over uj x using the 
identity 



duj- 



1 



= 2tt 



6(Re(A)) - 0(Re(B)) 



iuj - Aiuj - B B - A 

for Re(A) ^ or Re(B) ^ 0. Thus, 

ZS ^ 



(13) 



v x3,x'l v 4x',2x 



i(u)i - uj 3 ) + E Sx , nx , - E Si 
In a similar fashion, we have 

^ — V x4,x'l V 3x',2x 

Q(E Sx nJ ~ Q{E s .n x ,) 



(14) 



and 



i{Uli — LO4) + E Sx , nx , — E SxT , 
BCS = 2E^2i Q <6 77 ' 

^( E s x n x ) - 0(-E B , n ,) 
X = - — 1 

i(wi + lo 2 ) - E 3 , n , - E Sx 



(15) 



(16) 



We then approximate the denominators as, for example, 
- U3 3 ) + E Sx , nxl - E Sxnm w E Sx , nxl - E Sx7lx since Uuj c 
is much larger than the frequencies at which we probe 
the system. In principle, however, these effective inter- 
actions become retarded if we do not drop the frequency 
dependence. This approximation eventually yields 



70 T/7«'i7 / «T/i8V>i8'7 
° ~ V x3,x'l V 4x',2x 

w 6{E Sxnx ) - 9{E S , n ,) 



E m 



E, 



(17) 



ZS' 



rrjP^'arra'-y' ,M 

v x4,x'l v 3x',2x 



E 

Q{E Sx nx) ~ (^(Es^n^,) 

E ~ E 



(18) 



and 



BCS 



1 t / 77',/3q 1/ 0'q',77' 

2 2-^i xx', 21 v 43, xx' 

„Q{ E s x n x ) - 6{-E s , n ,) 



-E s 



Eg 



(19) 



We can extract some physics from these formulas with- 
out calculating any numbers. Each interaction matrix el- 
ement {y24x'i from the ZS diagram, for example) comes 
with a factor of e 2 /e£o according to Eq. fcfy. Therefore, 

once a factor of Huj c is pulled out of the denominator, 

(2) 

the ZS, ZS', and BCS contributions to U4 3 . 21 are propor- 
tional to K,(e 2 /e£o), as expected. In the lowest Landau 
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level (N — 0), the ZS and ZS' terms vanish since there 
cannot be any hole excitations in the internal legs (the 
x'j' leg) since there are no filled Landau levels below in 
which to have virtual hole excitations. Therefore, only 
the BCS diagram renormalizes the two-body interactions 
in the lowest Landau level-incidentally, this is not the 
case for graphene since there all Landau levels N = are 
filled in which it is possible to excite virtual holes. 
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FIG. 3: One of the nine three-body diagrams. 



(3) 

We now turn to the three-body interaction ^654.321- 
At order k there is a three-body term that is generated 
in our expansion of the action. Fig. [3] shows one of the 
nine diagrams that contribute to the three-body term. 
The expression fo r the sum of all nine diagrams can be 



order of Grassman variables, apart from keeping track of 
signs. However, when we pass from the functional inte- 
gral to the operator formalism, we must be careful about 
the order of these operators. The order of these operators 
is determined by a point-splitting regularization: 



du:e lu]0+ C t(u)c(w) = J(t)c{t) 
due lu0 ~ c\u)c{uj) = c(t)c\t) 



(23) 



In other words, the order of operators at the same time 
is determined by whether the contour must be closed in 
the upper- or lower-half-plane. In the case of our three- 
body inter acti on, which is actually a retarded interaction, 
as in Eq. (20 1, this is determined by the sign of E Stctna .. 
When we make the approximation of neglecting uj x in the 
denominator of Eq. ( 20 1 to pass to the low-energy limit 



in Eq. (21 1, we can drop the c^-dependence, except that 



we must remember the operator ordering that it imposes. 
Therefore, the three-body term that we have computed 
is not really 



E 



v 6x,21 v 54,x3 f f f 

—~ 4 c 5 c 4 c 3C 2 Cl 



(24) 



written compactl y * ' as 



but is, instead, 



E' E 



E E 

n x ,s x — m x —0 7— 4, eye. perm. 



v 6x,21 v 54,x3 



(20) 



At low- frequencies, this leads to a three-body term: 



(3) 
1 654:321 



- E E E 

« x .s x =0 m E =0 7=T>4- 



E 



eye. perm. 



v 6x,21 v b4,x3 



(21) 



where we sum over all cyclic permutations of the indices 
labeled with (1, 2, 3) and (4, 5, 6). The prime on the sum- 
mation over n x and s x indicates that we do not include 
cases where s x = and n x = N simultaneously. 



B. Normal-Ordering the Three-Body Interaction 

There is one important subtlety that we now discuss. 
When we compute the three-body term, just discussed, 
that is induced at order k, we evaluate diagrams that 
correspond to the operator expression: 



1 

2! 



c 6 c x C 2 C 1 C 5 C 4 c x c 3 + C 5 C 4 c 3 c x c x c 6 c 2 c l 



(22) 



Here, we have subsumed spin, Landau level, sub-band, 
and orbital indices into a single index for brevity. In the 
functional integral, we do not need to be careful about the 



frCt \,f}a-rrf}'*i' 
6x, 21 v hi,x3 



2 ^ 



JJJ 



\®(Es x ,n x )c\c\c\ C 3 C 2 C 1 + ^(~^s x ,n x ) c 6 C 2 C l C 5 C 4 C 3 
+6{ — Es x ,n x ) c 6 c 5 C 4 c 3 C 2 C 1 — ^s, ,n, ) C 5 C 4 c 3 C 6 C 2 C 1 



(25) 



After we normal-order, this can be re-written as: 



3! 



(3) t t t 

<54;321 c 6 c 5 c 4 c 3C 2 Cl 

v x3,x l v 4x ,2x ^ 

E 

f>A/3',7a T ">a'7,/3A 
V -A —1 v 3x ,2x 



(-^ Sx ,n x )c4C 3 C 2 Ci 



TVAa',A7-rV ( 9'a / ,A7 
J- V - ^ v xx ,2l v 43 



^(E StB , n j44c 2 ci . (26) 



2 ^ 



Therefore, the renormalization of the two-body interac- 
tion takes the form 



1 

2! 



<^43;21 = ZS + ZS' + BCS + ^U 4 3 J ;21 



1 



■„~.(2) 



(27) 



where ^^43-21 lii the sum of the last three terms in 



Eq. (26 1. Note that, in the N = Landau level, the 



first and second two-body terms in Eq. ( 26 1 vanish and 
only the third term contributes. 
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(2) 

There is a simple check that shows that <5w4 3 . 2 i is neces- 
sary. Suppose that we compute two-particle scattering to 
second-order in the Coulomb interaction. We can do this 
either in the full theory (with all Landau levels kept) or in 
the low-energy effective theory of the N th Landau level. 
The answers must be the same either way. With the 
Su 43 . 21 above, this is the case. There is a simple physical 
interpretation for the renormalization of the two-body 
term that results from normal-ordering the three-body 
term. We can think of this term as arising by taking 
a three-body vertex and connecting one of the incoming 
lines to one of the out-going lines in Fig. [3] If the two lines 
are connected to the same vertex (for example, 1 and 6 
m Fig. m then we have a tadpole diagram, which van- 
ishes. (Two such lines should not have been connected 
anyway since the operators at a single vertex are normal 
ordered.) If the two lines are connected to different ver- 
tices, then the resulting diagram is of the form of the ZS, 
ZS', or BCS diagrams (see Fig. [2]) , but one of the internal 
lines remains in the 7V th Landau level, rather than one of 
the higher or lower Landau levels that is being integrated 
out. In other words, this contribution accounts for pro- 
cesses such as the following: (i) two electrons in the N th 
Landau level interact, (ii) one of the is excited virtually 
to a higher Landau level while the other remains in the 
N th Landau level, and (iii) they interact again and the 
excited electron falls back to the N th Landau level. 

One might worry that the result of such a process will 
depend on the couplings within the ./V th Landau level and 
its filling factor. The two-body terms that we have gen- 
erated in this subsection derive from normal-ordering the 
three-body terms, so the dynamics within the iV th Lan- 
dau level does not enter into their computation. These 
dynamics (and the filling of the N th Landau level) en- 
ters only when we attempt to solve the resulting effective 
Hamiltonian for the N th Landau levels. 

Another way to think about this is to recall that we 
have not integrated out states within the N th Landau 
level. Within a Wilsonian renormalization group scheme, 
which is essentially what we have adopted, all internal 
lines in diagrams are at high energies and all external 
lines are at low energies because we integrate out N' ^ N 
Landau levels but do not do anything to states in the 
N th Landau level, which will be dealt with later (by, for 
instance, exact diagonalization). In our case, this means 
that all internal lines are in the N' ^ N Landau levels 
while all external lines are in the N th Landau level. In 
such an approach, a diagram with k external legs in which 
p internal lines are in the N' ^ N Landau levels while q 
are in the 7V th Landau level arises in the following way. 
When we integrate out the N' ^ N Landau levels, we 
will generate a vertex with k + 2q legs. One of the terms 
contributing to this vertex will have p internal lines. If 
we were to solve the problem in the N th Landau level 
perturbatively, then there would be a diagram in which 
we took this (A; + 2<7)-leg vertex and connected q incoming 
lines to q outgoing lines. Therefore, this physics is present 
in the (k + 2g)-leg vertex. In our case, this means that 



two-body processes in which the intermediate state has 
one electron in the N th Landau level and one in an N' ^ 
N Landau level are present in our three-body vertex. 
However, to correctly account for them, it is crucial to 
order the operators in the three-body vertex correctly. 

Note that in Fermi liquid theory for electrons with no 
external magnetic field (B = 0) we never have to con- 
sider the type of virtual process which gives this contri- 
bution because momentum conservation does not allow 
processes in which one electron is scattered to momenta 
far from the Fermi surface while the other electron stays 
near the Fermi surface. In low Landau levels, by con- 
trast, there is no conservation law precluding such pro- 
cesses and they are not only present but, in fact, give 
a substantial contribution to the renormalization of the 
two-body interaction. 



Effective Hamiltonian 



(2) 

Once we have these new interactions, u\ 3 . 21 an d 

(3) 

iig 54 . 321 , we can calculate the Haldane pseudopotentials 

and their three-body generalizations^. We use vffi s (N) 
to denote the bare two-body pseudopotential for two elec- 
trons residing in the N th Landau level with relative angu- 
lar momentum M and total spin S. The O(k) correction 
(including the contribution discussed in the previous sub- 

(2) 

section) to this pseudopotential is denoted by <W^ f g(iV). 
The three-body pseudopotential for three electrons with 
relative angular momentum M (defined as the degree of 
the relative wavefunction)^ and total spin S is denoted 
by V^ S (N): 



v ms( n ) = 53<Af5|m 3 a / J m4^ / )(m 1 a,m 2 )8|MS) 



x u 



(2) 
43;21 



(28) 



^2{MS\m 3 a', mi/3'} (rrna, m 2 p\MS) 



x 5u 



VmU N ) = ^(M5|m 47 ',m 5 /3',m 6 a') x 

(m 1 a,m 2 /3,m 3 -f\MS)u'^ 4 . 321 (30) 

where ^2 indicates a sum over all m, and primed spin 
variables. 

We will drop the spin indices from the pseudopoten- 

(2) 

tials and write (N) 



V§HN) = V£> s (N) 
where it will be clear from context what S is, i.e., for 
the two-body case with spin polarized (unpolarized) elec- 
trons S = 1(0) while for the spin-polarized three-body 
case S = 3/2. Our Hamiltonian is appropriate for all fill- 
ing factors. Therefore, we calculate even and odd M for 
the two-body terms in addition to M = 1, 2, 3, 5, 6, 7, 



0), 



8 



and 8 for the three-body terms. It should be clear what 
the total spin is for those particular M's. As shown in 
Section hill the numerical value for Vg (0), Vg (1) are 
very small, so we ignore M > 8. This allows us to avoid 
the potential complication of pseudopotential matrices, 
instead of pseudopotential numbers, that arise for total 
relative angular momentum M > 9 due to the multiple 
three-particle wavefunctions-^R 

Starting with Eq. ([2]), we have chosen to calculate the 
effective interactions and the resulting pseudopotentials 
in the planar geometry, rather than in the spherical ge- 
ometry. For the two-body pseudopotentials, we are only 
reporting the corrections due to Landau level and sub- 
band mixing and one must consider these corrections 
as corrections to the planar bare pseudopotentials. In 
order to calculate them in the spherical geometry, we 
would have to recalculate the corrections for each system 
size, which is extremely cumbersome and, presumably 
unnecessary (when we extrapolate to the thermodynamic 
limit). The spherical pseudopotentials approach the pla- 
nar pseudopotentials in the therm odyna mic limit of a 
large sphere and it has been arguecP^MJ t na t using the 
planar pseudopotentials in finite systems better approx- 
imate the infinite system. 

With these new effective pseudopotentials, we con- 
struct a Hamiltonian for the FQHE in the N th Landau 
level which includes Landau level and sub-band mixing 
effects perturbatively to first order in k, 



We also briefly comment about the projection operators: 
what is often not explicitly written is that the m-body 
projection operators Pij... m (M) have normalization set 
by the property that their eigenvalues must be or 1. 

Since every Landau level has the same number of 
states, we project the N th Landau level to the N = 
Landau level. The difference between the iV th Landau 
level and the zeroth is reflected in the values of the pseu- 
dopotentials. However, from the point of view of any cal- 
culations, our effective Hamiltonian, Eq. (31 ), is a Hamil- 



tonian for electrons in the lowest Landau level that sim- 
ulates the Hamiltonian for electrons in the N th Landau 
level. 



III. THE EFFECTIVE TWO- AND 
THREE-BODY PSEUDOPOTENTIALS FOR GaAs 

We now discuss the numerical values of the two- and 
three-body pseudopotentials in the presence of Landau- 
level and sub-band mixing appropriate to GaAs (i.e. for 
parabolic B = bands). For the two-body case, we 
calculate Eq. (27). The result is a function of single- 



H eB ( K ,d/£ ) = V^ 2) (d/£ ) + KVf>(d/£o) 



(2) 



+ KV 1 {3) (d/£ ) (31) 



particle angular momenta mi, m,2, JTI3, and 7714, denoted 
(2) . 

by 5u 43 . 2 it which depends on the Landau level index N 

of the partially-filled Landau level which we retain. All 

other Landau levels are either completely empty or com- 

(2) 

pletely full and are integrated out. Then, 5u 43 . 2 i is sub- 
stituted into Eq. (29), yielding the two-body pseudopo- 
tential correction 6Vffi(N). For the three-body term we 



where 



(Eq. (20)) which is a function of angu- 
and mi. The result, 



n 2 \d/(, ) = Y J v{\v l -v ] 



calculate Ug54-32i 

lar momenta mg, m.5, 7774,' 7713, 777,2, 
which depends on the index TV of the partially-filled Lan- 
dau level which we retain, is substituted into Eq. (30), 

yielding the three-body pseudopotential vjj (N) . 



= J2 V m ( N > d l 1 ^ ^' ( M ) ( 32 ) over Landau levels and sub- 



However, Eq. (27), and (20) all contain infinite sums 



A I 



i<3 



is the bare pseudopotentials with V(|rj — |) being 
the usual finite-thickness augmented Coulomb interac- 
tion (see Refs. 134133]) ; it is a function of the quantum- well 
thickness dfl§. In this expression, Pij(M) is an operator 
that projects the pair of electrons (i, j) onto a two-body 
state of relative angular momentum M. Similarly, 

Vf } (d/io) = E SV m (N, d/io) J2 h (M) (33) 



A I 



i<3 



and 



v} 3) (d/l )=Y,vP(N,d/e ) £ P ijk (M) (34) 

L i<j<k 

where Pijk(M) is a projection operator that projects onto 
triplets of electrons with relative angular momentum M. 
The subscripts and 1 in Eq. (31) indicate the bare 
interaction and the first order corrections, respectively. 



mnds and none of these sums 
can be done analytically. Therefore, they must be trun- 
cated and the resulting finite sums must be computed 
numerically. The infinite sum is the limit of the finite 
sums as the truncation is taken to infinity. We calcu- 
late our final results (5V$(N,d/t ) and vjp(N,d/£ )) 
as a function of the truncation and determine the con- 
vergence to the infinite sum by fitting to some chosen 
general common function. 

The two-body pseudopotential corrections and three- 
body pseudopotentials are given in, respectively, Tables 
[T] and |lT| for d/£ — 0, 1, 2, 3, and 4 and shown graphi- 
cally in Fig. [4] up to d/io = 8. The most salient feature 
is that these corrections are small compared to the bare 
pseudopotentials in the N — Landau level (except for 
M = 0) and the N = 1 Landau level (except for M < 3. 
This was already noted in Ref. [SHI for d/£o =0, however, 
without taking into account the modified two-body con- 
tribution. (Note that the numerical values which we find 
for d/£ — are different than in Ref. [591 due primar- 
ily to the correction to the two-body term coming from 
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FIG. 4: (Color Online) The Landau level mixing induced corrections to the two-body pseudopotentials for M = 1-6 in the lowest 

(N = 0, left figures) and second (N = 1, right figures) Landau levels. The top left plot in both panels shows all 5Vj^\N,d/lo) 
versus d/lo on the same plot for an easier comparison. Note that for the lowest Landau level (N = 0) the M — plot is on 
a different vertical scale than all the others since its value is so much larger. For the N = Landau level, the corrections are 
less than 10% (for all d/lo) of the bare pseudopotential value-discounting M = which is nearly 30% (for d/£o = 0) the bare 
pseudopotential. For the N = 1 the corrections are much larger. For M = the correction is nearly 60% and remains near 
50% for M < 3. Of course, for both Landau levels the corrections are mitigated by finite width (increasing d/lo). All energies 
are given in units of e 2 /clo- 



the correct normal-ordering of the three-body term.) By 
"small compared to the bare pseudopotential" , we mean 
that the coefficients (for M > in the N — Landau 
level and M < 3 for the N = 1 Landau level) are ss 10 
times smaller and for increasing angular momentum M 
approximately 50 times smaller. The smallness of these 
coefficients means that it is conceivable that our perturba- 
tive calculation is valid even for n beyond one. However, 
in order to say this with certainty, we would need to com- 
pute the order terms and show that the series 
converges. 

From Fig. |4j we can see the dependence of the two- 
body corrections to the bare pseudopotentials on d/lo. 
We observe in both the lowest and first excited Landau 
levels that generally the magnitude of the corrections to 
the two-body pseudopotentials decrease with increasing 
d/^o - remember that these numbers will be multiplied by 
K when constructing the full two-body pseudopotentials. 
The bare Coulomb interaction is itself also softened as 
d/£o increases and the fractional change in the bare pseu- 
dopotential and its two-body correction are similar. This 



behavior begins to break down very slightly in the N = 1 
Landau level at larger values of M where some mild non- 
monotonic behavior is observed, however, we emphasize 
that the relative value of the correction is quite small for 
higher M's and this non-monotonic behavior will most 
likely have little effect on the physics. The two-body 
corrections decrease with increasing d/£o and, eventually 
are expected to vanish, so in essence, the effect of finite- 
thickness is to mitigate the effect of Landau level mixing. 

In Fig. [5] we show the two-body corrections to the bare 
pseudopotentials as a function of relative angular mo- 
mentum M for d/l = 0, 1, 2, 3, and 4 for N = and 
N = 1. As in the case of the bare pseudopotentials we 
observe that the two-body corrections decrease with both 
increasing M and d/lo- However, an effect not seen in 
the bare pseudopotentials for N = 1 is observed in the 
two-body corrections; they are larger in magnitude than 
their N = counterparts. 

Since there is no three-body term at tree-level, the 
O(k) correction is the leading three-body term. Al- 
though the three-body interaction is numerically small, 
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TABLE I: We list the numerical values of 6V^\N, d/£ ) for 
d/£ = 0, 1, 2, 3, 4 and N = 0, 1. Values for d/£ = were 
given previously in Ref . [59] for M = 1 , 3 but due to the new 
two-body correction coming from properly normal-ordering 
the three-body interaction the values have been changed-^. 
All energies are given in units of e 2 /e£o- 







d/£ = 


1 


2 


3 


4 




(0,d/£ ) 


-0.3422 


-0.1963 


-0.1281 


-0.0900 


-0.0665 




(0,d/£ ) 


-0.0328 


-0.0300 


-0.0254 


-0.0211 


-0.0175 


5V} 2) 


(0, d/to) 


-0.0112 


-0.0108 


-0.0098 


-0.0088 


-0.0077 


svw 


(0,d/£ ) 


-0.0055 


-0.0054 


-0.0051 


-0.0047 


-0.0043 




(0,d/£ ) 


-0.0033 


-0.0032 


-0.0031 


-0.0029 


-0.0027 




(0,d/£ ) 


-0.0022 


-0.0022 


-0.0021 


-0.0020 


-0.0019 


sv™ 


(0,d/£ ) 


-0.0016 


-0.0015 


-0.0015 


-0.0014 


-0.0014 


sv 7 m 


(0,d/£ ) 


-0.0012 


-0.0011 


-0.0011 


-0.0011 


-0.0010 


svF> 


(0,d/£ ) 


-0.0009 


-0.0009 


-0.0009 


-0.0009 


-0.0008 


svf> 


(0,d/£ ) 


-0.0007 


-0.0007 


-0.0007 


-0.0007 


-0.0007 


5V (2) 


(l,d/£ ) 


-0.3816 


-0.3184 


-0.2696 


-0.2307 


-0.1997 




(l,d/£ ) 


-0.2143 


-0.2020 


-0.1815 


-0.1617 


-0.1442 


sv 2 {2) 


(l,d/t ) 


-0.1787 


-0.1456 


-0.1296 


-0.1172 


-0.1067 


5V^ 2) 


(l,d/£ ) 


-0.1039 


-0.0986 


-0.0927 


-0.0868 


-0.0812 


5Vi 2) 


(l,d/£ ) 


-0.0789 


-0.0729 


-0.0694 


-0.0664 


-0.0635 


SV^ 


(l,d/t ) 


-0.0353 


-0.0423 


-0.0458 


-0.0473 


-0.0476 


BVF 


(l,d/£ ) 


-0.0258 


-0.0296 


-0.0334 


-0.0361 


-0.0376 


sv 7 {2) 


(l,d/£ ) 


-0.0115 


-0.0181 


-0.0234 


-0.0272 


-0.0297 




(l,d/t ) 


-0.0073 


-0.0123 


-0.0175 


-0.0216 


-0.0245 


svf> 


(l,d/£ ) 


-0.0023 


-0.0079 


-0.0133 


-0.0176 


-0.0209 



as may be seen in Table [II] and Figs. [6] and [7| it is the 

leading term which breaks particle-hole symmetry. The 
three-body term is especially interesting for the FQHE at 
v = 5/2. In its absence, particle- hole symmetry is a sym- 
metry of the Hamiltonian and, consequently, the most 
promising candidates to describe the the 5/2 FQHE, the 
Moore-Read Pfaffian state and the anti-Pfaffian state, 
would be exactly degenerate. As is well knowrP2E3 ; the 
Moore-Read Pfaffian is the exact zero-energy solution to 
a repulsive (positive) three-body Hamiltonian that only 
contains an M = 3 term. The anti-Pfaffian, of course, is 
the exact zero-energy solution to the Hamiltonian that 
is the particle-hole conjugate of the one that yields the 
Moore-Read Pfaffian. Under particle-hole conjugation, 
the positive M = 3 three-body term picks up a minus 
sign and becomes attractive. Additionally, however, a 
two-body term is also generated under particle-hole con- 
jugation of the Moore-Read Pfaffian Hamiltonian. The 
M = 3 three-body pseudopotential due to Landau level 
mixing is negative in both the N = and N = 1 Landau 
levels, a sign which, naively, favors the anti-Pfaffian state 
over the Moore-Read Pfaffian. However, with the M — 5, 
6, 7, and 8 three-body terms being non-zero, one cannot 
predict the physics from the pseudopotentials alone. 

To illustrate the d/£o dependence of the three-body 



TABLE II: We list the numerical values of the three-body 
pseudopotentials V^'(N, d/£o) as a function of M for d/£o = 
0, 1, 2, 3, and 4 and N = 0, 1. The values for d/l = were 
given previously in Ref. 1591 where again a few values have been 
changed very slightly due to a more careful extrapolation. All 
energies are given in units of e 2 /e£o- 







d/£ = 


1 


2 


3 


4 


v} 3) 


(0,d/£ ) 


-0.0345 


-0.0400 


-0.0364 


-0.0312 


-0.026k3 




(0,d/£ ) 


-0.0540 


-0.0410 


-0.0324 


-0.0262 


-0.0216 


v 3 (3) 


(0,d/£ ) 


-0.0181 


-0.0173 


-0.0156 


-0.0138 


-0.0120 




(0,d/£ ) 


0.0033 


0.0026 


0.0015 


0.0006 


-0.0000 


T/( 3 ) 


(0,d/£ ) 


-0.0107 


-0.0102 


-0.0093 


-0.0083 


-0.0073 


v 7 (3) 


(0,d/£ ) 


0.0059 


0.0054 


0.0043 


0.0033 


0.0025 


y(3) 
V S 


(0,d/£ ) 


-0.0047 


-0.0045 


-0.0041 


-0.0037 


-0.0033 


vf> 


(I, d/to) 


-0.0319 


-0.0280 


-0.0232 


-0.0192 


-0.0160 


vf> 


(l,d/£ ) 


-0.0305 


-0.0223 


-0.0174 


-0.0142 


-0.0117 




(l,d/£ ) 


-0.0147 


-0.0136 


-0.0118 


-0.0101 


-0.0087 


vf> 


(l,d/£ ) 


-0.0054 


-0.0051 


-0.0047 


-0.0042 


-0.0038 


v 6 (3) 


{I, d/to) 


-0.0099 


-0.0093 


-0.0082 


-0.0071 


-0.0061 


vf' 


(l,d/£ ) 


0.0005 


0.0002 


-0.0001 


-0.0004 


-0.0006 


y(3) 


{l,d/£ ) 


-0.0009 


-0.0013 


-0.0016 


-0.0017 


-0.0018 




FIG. 5: (Color Online) SVffi (AT, d/to) as a function of M 
for d/£o = 0, 1, 2, 3 and 4. The left panel corresponds to the 
lowest Landau level (N — 0) while the right panel corresponds 
to the first excited Landau level (N — 1). As in the case of 
the bare pseudopotentials, the corrections to the two-body 
pseudopotentials decrease with increasing M. All energies 
are given in units of e 2 /e£o- 



pseudopotentials, we show V^(N, d/£o) as a function of 
d/£o and M in, respectively, the N — and N = 1 Lan- 
dau levels in Figs. [6] and [7] Other than the M = 1, 
N = case (discussed below), finite width decreases 

V^iN.d/^) monotonically. Fig. M shows V^\N,d/£ Q ) 
as a function of M and we see that for M > 5 the ab- 
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FIG. 6: (Color Online) The three-body pseudopotentials 

Vm (N, d/lo) as a function of d/lo for the lowest Landau level 
(left panel) and the first excited Landau level (right panel). 
Note that these are numerically small compared to the bare 
two-body pseudopotentials. All energies are given in units of 
e 2 /el . 




FIG. 7: (Color Online) V^'(N,d/l ) as a function of M for 
d/lo = 0, 0.5, 1-6 for the lowest Landau level (left panel) 
and the first excited Landau level (right panel). The magni- 
tude decreases with M except for the non-trivial oscillatory 
behavior for M > 5. All energies are given in units of e 2 /e£o- 



solute values oscillate: \Vf> {N,d/£ )\ < \V^> {N,d/£ )\, 
\V 7 (3) (N,d/£ )\ < \V^(N,d/£ )\, etc. 

We note that there is some non-trivial behavior in the 
M = 1 and 2 three-body pseudopotentials for small val- 
ues of d/£o for N = 0, 1. The M = 1 pseudopotential 
in the lowest Landau level is non-monotonic and starts 
out smaller than the M = 2 pseudopotential. As d/lo 
increases, the two cross and M = 1 is larger than M = 2 



for large d/£o. In the N — 1 Landau level, we see the 
same qualitative behavior except that M = 1 and 2 do 
not change signs. 

Finally we examine the overall structure of the effective 
interaction by showing vffl (N, d/lo) for N = 0, 1 (top 
and bottom panels, respectively) and p — 2, 3 in Fig. [8] 
as a function of M for a range of ft. The two-(three-)body 
pseudopotentials are color-coded with "hot" ( "cold" ) col- 
ors going from red to yellow (green to blue) as k is 
changed from to 3. The numerical values of the cor- 
rections to the bare two-body pseudopotentials and the 
three-body pseudopotentials are also small compared to 
the bare two-body pseudopotentials. However, as noted 
in Sec. |TJ small quantitative changes to the pseudopo- 
tentials can have large qualitative effects on the possible 
ground states of the effective FQHE Hamiltonian. As 
noted above, naively one might expect that the negative 
value of the M = 3 three-body term might favor the 
formation of the anti-Pfaffian, but for values of K > 1 
there is a fair amount of renormalization of the two-body 
pseudopotentials. It is known from previous worliMMH 
that finite-thickness stabilizes the Moore-Read Pfaffian 
at v = 5/2 but this calculation was only done for k = 
and could not distinguish between the Moore- Read Pfaf- 
fian and anti-Pfaffian. From examining the Landau level 
and sub-band mixing induced corrections to the Haldane 
pseudopotentials there is potentially a range of d/lo and 
k that would favor or disfavor the Moore-Read Pfaffian 
and anti-Pfaffian. In fact, a cursory look at the two-body 
pseudopotentials lead one to suspect that Landau level 
mixing is relatively unimportant for spin-polarized elec- 
trons in the N = Landau level, while evidently, Lan- 
dau level mixing is much more important in the N = 1 
Landau level. However, we emphasize that the effect of 
the three-body terms on the physics of the FQHE in ci- 
ther the N = or 1 Landau levels is non-trivial and not 
known. We are currently undertaking a thorough numer- 
ical study to determine the full quantum phase diagram 
of this new Hamiltonian for GaAs^H 

We now turn our attention to the study of an effective 
Hamiltonian for graphene. 



IV. EFFECTIVE HAMILTONIAN FOR 
GRAPHENE 

The formalism for including Landau level mixing in 
graphene (linear dispersion at B — 0) is similar to that 
of GaAs, so we only highlight the differences, such as 
the absence of finite thickness. The Landau level prob- 
lem has been investigated by many authoriP^^ anc j we 
only briefly repeat those results that are relevant to the 
present discussion. 

The action (Eq. (1)) for graphene is of the same form 
as for GaAs with three differences. The first, and most 
important, is that the kinetic energy of an electron in 
the A^ th Landau level is E n = sgn(n) v/ '2|n|?iw p / '£ a where 
vf is the electron Fermi velocity (~ 10 6 m/s). Unlike 
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FIG. 8: (Color Online) The two-body and three-body pseudopotentials, (TV, d/lo) and (TV, d/lo), as functions of M for 

k = 0-3 for TV = (top panels), TV = 1 (bottom panels) for d/l = 0(left), 2 (middle), and 3(right). Obviously, V$ (TV, d/£ ) is 
a sum of the bare two-body term plus k times the two two-body corrections. The two-body pseudopotentials are color-coded 
such that k = corresponds to red and k — 3 corresponds to yellow, i.e., the color-coding is "hot". The three-body terms are 
color-coded "cold", i.e., k = is green and the color changes continuously to blue for k — 3. Generally, increasing d/lo and 
k reduces the two-body pseudopotentials compared to the bare values (« = 0). The three-body pseudopotentials are directly 
proportional to k so increase in magnitude for increasing k, however, increasing d/lo mitigates this effect by decreasing the 
coefficients themselves, as expected. Please note that any two-body corrections for M > 10 are expected to be small and 
unimportant. From inspection of this figure one would expect that Landau level mixing is likely to be unimportant for the 
FQHE in the TV = Landau level for spin-polarized electrons. For the TV = 1 Landau level, however, it is not clear what effect 
these changes to the pseudopotentials will have on the physics and the answer will have to await future exact diagonalization 
studies. Of course, even in the TV = Landau level the effect of the three-body terms cold produce non-trivial effects. All 
energies are given in units of e 2 /elo- 



in GaAs, this kinetic energy can be positive or negative 
(electrons or holes) and increases as the square root of the 
absolute value of the Landau level index n. Hence, the 
Landau levels in graphene are not equally spaced; the 
spacing between successive Landau levels decreases (as 
the inverse square root) as the Landau level index is in- 
creased. The second important difference is that there is 



no sub-band index "s" since the thickness of a graphene 
monolayer is atomically thin. The third and final dif- 
ference lies in the interaction matrix elements. We will 
use the notation for the matrix elements of the 

Coulomb interaction for electrons in graphene and we 
will continue to use 143,21 for the corresponding matrix 
elements in GaAs. 



In GaAs, the electron-electron interaction matrix element 143,21 is given in Eq. (3). In graphene, Eq. (3) is modified 
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/d 2 k 2 \/2^~' i=1 Sni '° 

~2^V(k)e k - G mim2 (—k)G m3mi (k) 

x{sgn(n4)sgn(n 2 )G| n4 |_i,| n2 |_i(-fc) + G|„ 4 ||„ 2 |(-fc)} 

x{sgn(n 3 )sgn(ni)G| n3 |_ 1) | ni |_ 1 (fc) + G\ nz \\ ni \(k)} , (35) 

where we set sgn(O) = in the above equation. In this equation, V(k) = e 2 / (2ire£ok) is the pure (i.e. zero-thickness) 
Coulomb interaction. Note that the definition of Vf^^ 1 follows from the the single-particle energy eigenstates r]f^ ph , 
which are related to the single-particle eigenstates r\ nm of GaAs according to: 



l^grapM 
I Inra I 



l 



I -sgn(n>'|r?| n |_ 1)TO ) \ ^ 
V%J \ \n\n\m) j 



Hence, if n» = for i = 1, 2, 3, 4, then V^2i = V43.21. Thus, the bare pseudopotentials for graphene in the N = 
lowest Landau level are identical to those in zero-thickness GaAs. In the ./V th Landau level, the pseudopotenials in 

graphene are a mixture of the N th and (N — l) th zero-thickness GaAs Landau levels. 

I 



Furthermore, the Landau level index N in graphene is 
not bounded from below by zero but instead spans — 00 
to +00. Hence, the sum in Eq. (11) is modified in the 
case of graphene to 



E = E E' E 



(37) 



>7>7'=4-,T 



Again the prime on the sums over n x and n x i indicates 
that we do not include n x = N or n x i — N in the sums. 

Another difference between graphene and GaAs is the 
Landau level mixing parameters. The Landau level mix- 



ing parameter for graphene is k 



2 /ee )/(hv F /£ ) 



e /(evph) = a g where a g is the graphene fine-structure 
constant. Notably, k is independent of magnetic field 
strength B. This is in stark contrast to GaAs where 
k oc l/\/~B, so that Landau level mixing can be ignored 
for sufficiently large B. The Landau level mixing pa- 
rameter in graphene k is a constant that can only be 
adjusted by manipulating the dielectric constant e, the 
Fermi velocity vp, or Ti\ In freestanding (i.e., suspended) 
graphene k ~ 2.2^^ while for graphene on substrates it 
is k fa 0.96 for a Si0 2 substrateP ( e pa 2.3) and k ps 0.44 
for a boron nitride substrate^ (e ps 5). Thus, by ad- 
justing the dielectric constant of the substrate from high 
dielectric to vacuum, one can only tune the Landau level 
mixing parameter between < k < 2.2 but it cannot be 
continuously tuned within a given device by tuning an 
experimental "knob" such as the magnetic field strength. 

We also emphasize that since the spacing between Lan- 
dau levels varies, unlike in GaAs where it is constant, and 
we expect the effect of Landau level mixing will depend 
on the Landau level in a manner that varies approxi- 
mately as cx k\J N. 



V. THE EFFECTIVE TWO- AND THREE-BODY 
PSEUDOPOTENTIALS FOR GRAPHENE 

We first investigate the three-body pseudopotentials in 
graphene since the discussion will help facilitate the later 
discussion of the two-body terms. 

The most striking feature is that the three-body pseu- 
dopotentials vanish for N = 0. Recall the expression for 
the three-body term generated in the perturbative ex- 
pansion of the action for GaAs (Eq. p0|)). For clarity, 
we write only the sum over Landau level index n x (and 
suppress the sums over m x , spin, and cyclic permutations 
since they are not relevant for this argument) 



E 

— 00,71^^0 



v 6x,21 v 54,x3 



^_ 
hvp 



E 



,->q'A,/3qtV^'7',A7 
v 6x,21 v 54,x3 

sgn(n x ) v /2|n 2; | 



E 



v 6x,21 v 54,x3 



n x — — oo 



E 



v 6x,21 v 54,x3 



E 



T>a'A,^af>/9'7',A7 
v 6x,21 v 54,x3 



f/Qt'A,/9af/ J 8 / 7',A7 
v 6x,21 v 54,x3 



= 

where Vf 3 %{^ a 

= E n:D 

sgn(n)^/2|n| 



V 3i2 i5 aP 'Sf Sa ' and 



V i3 , 21 S aa 8ft 

- fi = (hvp/£ ) (sgn(n x ) y/2 \n x | - 
The canceling of the two terms in the 
second to last line above can be seen by going back to 
the definition of V43,2i- Thus, the three-body terms ex- 
actly vanish in the lowest Landau level (N — 0). In the 
N = Landau level in graphene, particle-hole symmetry 
is an exact symmetry (in the absence of disorder) be- 
cause there are as many Landau levels below N = as 
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TABLE III: We list the numerical values of the three-body 
pseudopotentials vj^ (N) as a function of M for graphene. 
All energies are given in units of e 2 /e£o- 



TABLE IV: We list the numerical values of SV^(N) for 
N = 0, 1, and 2. The two-body correction consists of two con- 
tributions: one is a direction calculation of the ZS, ZS', and 
BCS diagrams and the other ari ses from the normal ordering 
of the three-body term (Eq. (271, please see the discussion in 
the text above. In the N = Landau level, these corrections 
are similar in magnitude to the results for GaAs. However, 
for N — 1 and 2 the corrections are much larger and, in fact, 
it is expected that effective Landau level mixing parameter k 
is really oc kvN, see discussion in text. We use the N — 
values as the first term contribution for all N since this term 
is so much smaller in magnitude than the second term and 
expected to make little difference to the final values or the 
physics. All energies are given in units of e 2 /e£o. 



N = 





1 


2 


3 


4 


V} 3) (N) 





-0.1237 


-0.0600 


-0.0377 


-0.0272 


% (3) (N) 





-0.0856 


-0.0556 


-0.0363 


-0.0271 


V 3 (3) (N) 





-0.0537 


-0.0413 


-0.0266 


-0.0194 


V 5 (3) (N) 





0.0135 


-0.0225 


-0.0221 


-0.0165 


K (3) (N) 





0.0313 


-0.0262 


-0.0225 


-0.0164 


V 7 (3) (N) 





0.0205 


-0.0105 


-0.0133 


-0.0141 


V S (3) (N) 





-0.0123 


0.0027 


-0.0148 


-0.0147 



above. Consequently, Landau-level mixing cannot gener- 
ate a 3-body term because it would violate particle-hole 
symmetry. 

Outside of the lowest Landau level (N ^ 0) the three- 
body pseudopotentials are non-zero. They are equal in 
magnitude and opposite in sign in the N th Landau level 
compared to the — N th Landau level (see Fig. [9]). This, 
again, follows from particle-hole symmetry. Furthermore, 
they are, in general, negative for N > and, therefore, 
positive for N < 0. 

In fact, the Hamiltonian for graphene is invariant un- 
der the combined action of sublattice symmetry, time- 
reversal, and charge-conjugation. The sublattice symme- 
try takes c — > — c on one sublattice, thereby inverting the 
kinetic energy in the absence of a magnetic field. Charge- 
conjugation transforms particles into holes, thereby in- 
verting the kinetic energy in the absence of a magnetic 
field but leaving the Coulomb interaction energy un- 
changed (up to a shift of the chemical potential). Time- 
reversal inverts the direction of the magnetic field. There- 
fore, the combination of all three transformations leaves 
both the kinetic and Coulomb energies unchanged. The 
combination transforms electrons in the N th Landau level 
into holes in the — 7V th Landau level. Therefore, we only 
give N > results. The N < values may be obtained 
by symmetry, as described above. 

Quantitatively, the three-body terms in graphene are 
similar in magnitude to those in GaAs, small compared 
to the bare two-body pseudopotentials, and decreasing in 
magnitude with increasing M. For \N\ = 1 there is non- 
trivial, non-monotonic M-dependence, where the pseu- 
dopotentials for M — 5, 6, and 7 are opposite in sign to 
those for M = 1, 2, 3, and 8. Table [TTT] lists the numerical 



values of vffl(N) for \N\ = - 4 and M = 1 - 3, and 
5-8. 

Turning now to the 2-body pseudopotentials, we note 
first that because the energy increases with Landau level 
index n as y/n, the contribution from higher levels will 
be less strongly suppressed. In fact, unlike in the case of 
GaAs, it is extremely helpful to use the constraint follow- 
ing from angular momentum conservation to eliminate 
one of the Landau level index summations in Eqs. (14)- 
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9: (Color Online) (A) as a function of M for N 



(16 1, thereby making the convergence of the sum more 



0-6. Note that for N = all three-body pseudopotentials 
exactly vanish. All energies are given in units of e 2 /e£o- 



apparent. 

The contribution to the two-body term coming from 
the normal ordering is usually an order of magnitude 
larger than the term coming from the ZS, ZS', and BCS 
diagrams alone. Because of this fact, and the fact that 
finding the convergence of the smaller two-body terms 
is laborious, we use the N — values of the two-body 
terms for all Landau levels as an upper limit and good 
approximation. Table IV lists the numerical values of 
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FIG. 10: (Color Online) The two-body pseudopotential cor- 
rections, SVjjj {N) for N — 0, 1, and 2 for graphene as a func- 
tion of M. It is clear from this figure that the pseudopotential 
corrections increase for increasing Landau level N reflecting 
the conjecture that the effective Landau level mixing param- 
eter is actually oc Ry~N and not simply k. All energies are 
given in units of e 2 /elo- 



V { m\n). 

In addition, in Fig. 10 we plot the Landau level mixing 
induced corrections to the two-body pseudopotentials as 
a function of relative angular momentum M for N = 0, 1, 
and 2. The two-body corrections are universally negative 
and for small M can be approximately half the value of 
the bare pseudopotential. Interestingly, the corrections 
become larger in higher Landau levels owing to the fact 
that the effective Landau level mixing parameter is oc 
K\/N. 



In Fig. 11 we show the two- and three-body pseudopo- 
tentials as a function of both k and M. For N = 0, 
Landau level mixing has little effect on the pseudopoten- 
tials, i.e., the corrections to the two-body term are rela- 
tively small and there is no three-body term-this might 
help explain why the FQHE is observed in the N = 
Landau level in graphene, when the system is spin and 
valley polarized, even though Landau level mixing cannot 
naively be ignored due to a large value of k. This further 
underscores the fact that the FQHE in graphene in the 
N = is observed to be, and is expected to be, nearly 
identical to that of GaAs (other than interesting compli- 
cations that arise when the spin and valley degeneracies 
remain). However, for N ^ 0, there are strong correc- 
tions in the form of negative two-body corrections and 
emergent three-body terms. For moderate values of k it 
is not at all clear that graphene would even exhibit the 
FQHE according to our calculations-the two-body terms 
can be of the same order of magnitude as the three-body 
terms. We also emphasize that with such strong Landau 
level mixing corrections to first order in k it is possible 
that higher order corrections will be significant, however, 
it is likely that these corrections (which we are ignoring) 



will make the situation worse for the FQHE, not better. 

Of course, it remains to be seen theoretically what ef- 
fect the three-body terms have on the FQHE. Exact diag- 
onalization studies using the above effective Hamiltonian 
are likely to shed light on this problerrPSI. 



VI. CONCLUSION 

Quantum electrodynamics is a paradigm of a pertur- 
bative theory: the dimensionless expansion parameter 
a ss 1/137 is small, so effects that occur at each order in 
perturbation theory are successively less important than 
those that occurred at the previous order. However, in 
condensed matter physics, we have grown accustomed to 
perturbation theory in dimensionless parameters that are 
not particularly small. Landau parameters are typically 
0(1), yet Fermi liquid theory remains valid for 3 He and, 
often, for electrons in metals. The ratio of the Coulomb 
energy to the kinetic energy in a metal, r s , can be moder- 
ately large, without the metallic state being destabilized 
(in favor of, say, a Wigner crystal) . But we should guard 
against the possibility that we have been lulled into a 
false sense of security. There is no quantum Hall exper- 
iment in which Landau-level mixing is obviously negligi- 
ble a priori. In the original observation of the fractional 
quantum Hall effect at v — 1/3, the Landau level mixing 
parameter is k = 0.65. This might be small enough that 
one can reach this value via an expansion in powers of k 
about k = or it might not, but it is certainly not a very 
small number such as 1/137 which one a priori expects 
to be negligible. In the N = 1 Landau level, 1 < k < 2.5, 
so it is even less clear that Landau level mixing is small. 
In the case of graphene, 0.44 < k < 2.2, so the situation 
is no better. 

The coefficient of O(k) corrections to two- and three- 
body pseudopotentials due to Landau level mixing are 
almost always numerically small. Consequently, these 
corrections are small even at K = 1. (Most of the cor- 
rections are small even at n = 10.) Since we have not 
computed the order k 2 corrections, this does not prove 
that experiments are effectively in a small n regime. How- 
ever, it does, at least, raise this possibility - or, in other 
words, the possibility that the true expansion parameter 
is - k/10. 

Of course, for the FQHE in the N ^ Landau level of 
graphene, the corrections are not necessarily small and it 
is very possible that a perturbative treatment of Landau 
level mixing is inappropriate. 

Strictly speaking, our computation yields an effective 
action, rather than an effective Hamiltonian - as is al- 
ways the case when integrating out high-energy states 
to produce a low-energy effective theory. This effective 
action has retarded interactions but, at lowest-order in 
k, these retarded interactions can be neglected. More- 
over, our calculation assumed fixed chemical potential. 
At lowest-order in k, we can simply take this to be the 
Landau level energy hu) c (N+ 1) or sgn(N) y/2\N\hvF / &o, 
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FIG. 11: (Color Online) Both the two-body and three-body pseudopotentials as functions of M for k — — 3 for graphene, in 
the N = 0, 1, and 2 Landau levels. For the lowest Landau level (N — 0) there is no three-body term (it exactly vanishes) and 
the two-body corrections are modest. Hence, one would expect, from merely studying the pseudopotentials, that the FQHE 
in N — Landau level for graphene would be nearly identical to that of GaAs. However, for N 7^ 0, we observe significant 
corrections to the two-body pseudopotentials and significant three-body terms. It is not at all clear what these effects will have 
on the FQHE. Naively, one is tempted to suggest that the FQHE will not be observed in N 7^ Landau levels in graphene, or 
if it is observed, it will be of an exotic variety. All energies are given in units of e 2 /ela. 



regardless of filling fraction within the N Landau level. 
However, in going to higher-order in re, we would have 
to tune the chemical potential to determine the filling 
fraction. For these two reasons, it will be difficult to de- 
termine the ground state of the system to order re 2 . But 
it should still be possible to derive an effective action (al- 
beit with retarded interactions) from which it would be 
possible to determine if order re 2 effects are also small 
even at re = 1. 

We note that, for large numbers of electrons, fixing the 
chemical potential will be the same as fixing the electron 
number (or, equivalently, the filling fraction). However, 
for the small numbers of electrons considered in numeri- 
cal exact diagonalization studies, fixed chemical potential 
and fixed electron number are not the same. For fixed 
chemical potential, An cx ^/n (where n is the electron 
number) which can be large. Thus, our effective Hamilto- 
nian is only, strictly speaking, correct for systems with a 
large number of electrons, not a small number. However, 
if the goal of studying small systems is to extrapolate 
their properties to larger ones to determine the solution 
in the thermodynamic limit with fixed filling fraction, 
then our effective Hamiltonian can be used. 

Even if it is quantitatively small, Landau-level mix- 
ing could, nevertheless, have significant qualitative ef- 
fects. It is the leading effect that breaks particle-hole 
symmetry within a Landau level (except for the N — 
Landau level in graphene, where this is an exact sym- 
metry). This has important effects at v = 5/2 in GaAs, 
where it could break the degeneracy between the Moore- 
Read Pfaffian state and the anti-Pfaffian state, as we 
discuss elsewhere^. It can also have important effects 
at v = 7/3,8/3 and 12/5, 13/5-with particle symmetry 



broken, these states may not be particle-hole conjugates 
of each other after all but, instead, have very different 
characters. Indeed, even within the lowest Landau level, 
there are obvious differences between particle-hole conju- 
gate plateaus, such as v — 1/5 and v = 4/5. Presumably 
the underlying cause is Landau level mixing. 

Finally, we note that, although the corrections to the 
effective Hamiltonian due to Landau level mixing are 
small compared to the bare terms, they are not small 
compared to the energy gaps of quantum Hall states. So 
it is possible, in principle, for three-body terms to have 
a non-negligible effect on the energy gaps of some frac- 
tional quantum Hall states while, at the same time, order 
re 2 terms can be neglected. In such a case, our order re ef- 
fective Hamiltonian would be the most promising starting 
point for an attempt to make a direct quantitative com- 
parison between experiments and numerical simulation 
of quantum Hall systems. 
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Note: during the final preparation of this manuscript 
we received a first draft of a paper by Sodemann and 
MacDonalcP containing some results on the effect of 
Landau level mixing on the FQHE in GaAs in the zero- 
thickness limit. They reproduced the d/£o — limit of 
our three-body interaction (which, in turn, is equal to the 
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three-body interaction of Bishara and Nayal<P^) and re- 
produced the d/£o — limit of our two-body interaction 
in the N = Landau level. However our two-body inter- 
action in the N = 1 Landau level diverges from theirs. In 



contrast, the <1/£q — limit of our two- and three-body 
interactions in both the N = and N = 1 Landau lev- 
els are within a few percent of the results of Rezayi and 
SimorPl 
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